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Abstract 

We find that for small vacuum mixing angle 9 and low energies (s <C M|) the width 
of matter, di/ 2 , needed to have conversion probability P > 1/2 should be larger than 
dmin = vr/(2^/2G F tan20): d 1/2 > d 

m in- Here Gp is the Fermi constant, s is the 
total energy squared in the center of mass and Mz is the mass of the Z boson. The 
absolute minimum — d m in is realized for oscillations in a uniform medium with 
resonance density. For all the other density distributions (monotonically varying 
density, castle wall profile, etc.) the required width dx/2 1S larger than d m i n . The 
width d m i n depends on s, and for Z-resonance channels at s ~ M§ we get that 
dmin(s) is 20 times smaller than the low energy value. We apply the minimum 
width condition, d > dmin, to high energy neutrinos in matter as well as in neutrino 
background. Using this condition, we conclude that the matter effect is negligible 
for neutrinos propagating in AGN and GRBs environments. Significant conversion 
can be expected for neutrinos crossing dark matter halos of clusters of galaxies and 
for neutrinos produced by cosmologically distant sources and propagating in the 
universe. 
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1 Introduction 



Since the paper by Wolfenstein the neutrino transformations in matter became one 
of the most important phenomena in neutrino physics. Neutrinos propagating in matter 
undergo coherent forward scattering (refraction) described at low energies by the potential 

V = V2G F n , (1) 

where Gp is the Fermi constant, and n is a function of the density and chemical compo- 
sition of the medium. For the case of u e — and u e — v T conversion in matter n coincides 
with the electron number density, n e . 

Refraction can lead to an enhancement of oscillations in media with constant density, 
and to resonant conversion in the varying density case (MSW effect) 0,0]. For periodic, 
or quasi-periodic density profiles, various parametric effects can occur [|]-|6|]. 

The MSW effect has been applied to solar neutrinos M, and to neutrinos from su- 
per novae 0. Oscillat ions of neutrinos of various origins (solar, atmospheric, supernovae 
neutrinos, etc.) in the matter of the Earth have been extensively studied. Apart from 
resonance enhancement of oscillations, parametric effects are expected for neutrinos cross- 



ing both the mantle and the core of the Earth [pHTljl. The oscillations and conversion of 
active neutrinos into a sterile species can be important in the Early Universe [O. Re- 
cently, matter effects on high energy neutrino fluxes from Active Galactic Nuclei (AGN) 



and Gamma Ray Bursters (GRBs) have been estimated [13]. Propagation of ultra-high 
energy neutrinos in halos of galaxies has been considered Jl4| . 

It is intuitively clear that to have a significant matter effect a sufficiently large amount 
of matter is needed. Let us define the width of the medium as the integrated density along 
the path travelled by the neutrino in the matter: 



d = J n e (L)dL . (2) 



This quantity is frequently named "column density" in astrophysical context. We will 
show that there exists a minimum value d m i n for the width below which it is not possible 
to have significant neutrino conversion. This lower bound is independent of the density 
profile and of the neutrino energy and mass. That allows us to make conclusions on 
the relevance of matter effect in various situations without knowledge of the density 
distribution. 

The paper is organized as follows: in section 2 we derive the minimum width con- 
dition for the conversion in matter between two active neutrino flavours, and check it 
for different density profiles. In section 3 we discuss the generalizations of the condition 
to the active-sterile case and to conversion induced by flavour changing neutrino-matter 
interactions. We also study the matter effect in the small width limits. Section 4 presents 
a study of the minimum width condition for high energy neutrinos both in matter and 
in neutrino background. Section 5 is devoted to applications of our results to neutrino 
propagation in AGN and GRBs environments, in dark matter halos and in the Early 
Universe. Conclusions and discussion follow in section 6. 
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2 The minimum width condition 



In this section we consider various mechanisms of matter enhancement of neutrino flavour 
conversion. For each of them we work out the minimum width of the medium needed to 
have significant conversion probability, showing that a lower bound for the width exists 
and is realized in the case of uniform medium with resonance density. 



2.1 The absolute minimum width 

Let us consider a system of two mixed flavour states^, v e and [y T \ characterized by 
vacuum mixing angle 9 and mass squared difference Am 2 . In a uniform medium the states 
oscillate and the transition probability, as a function of the distance L, is given by: 

P Ue ^(L) = sin 2 26 m sin 2 , (3) 

where 8 m and l m are the mixing angle and the oscillation length in the medium: 

. sin 20 

Sm m ~ [(2EV/Am? - cos2#) 2 + sin 2 26] 1 / 2 

m [{2EV/Am 2 - cos2fl) 2 + sin 2 20] " 1 ; 

Here I = AtiE/Am 2 is the vacuum oscillation length and E is the neutrino energy. 

We assume that the vacuum mixing is small, so that vacuum oscillations effects are 
negligible (P^^ v 1) and a strong transition in medium, i.e. P v ^ Vii = 0(1), is 
essentially due to matter effect. For definitness, we choose the condition of significant 
conversion to be 

P n >\- (5) 
Let us consider a uniform medium with resonance density || : 

Am 2 

res = COS 29 . (6) 

2V2EG F 

In this case the oscillation amplitude is sm28 m = 1, and the oscillation length equals 

jres = 1 = ^E 

sin 29 Am 2 sin 26 ' 1 ' 

According to eq. (§) the condition (|^) starts to be satisfied for L = l m /A, and the 
corresponding width is: 

A \ res ires (q\ 

^min 4 ' 



1 The arguments remain the same for three neutrinos. 
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Inserting the expressions of n r e es and l res given in @ and (0), we get: 

d ■ = - = (9) 

mm 2 v / 2G F tan2# tan2# ' v ' 

where 

d n ~ 111 (10) 
2 \p2G p Gf 

We will call do the refraction width. Numerically, 

d = 2.45 • 10 32 cm- 2 = 4.08 ■ 10 8 A cm" 2 , (11) 

where A = 6 • 10 23 is the Avogadro number] 2 ]. 

The widths d m i n and do have a simple physical interpretation. The refraction width 
d Q is a universal quantity: it is determined only by the Fermi coupling constant, and does 
not depend on the neutrino parameters at all. Using the definition of refraction length 

2tt _ 2tt 

/o = V ~ 7^g~ f (12) 

we can write: 

(13) 

n e 4 

It appears that c?o corresponds to the distance at which the matter-induced phase differ- 
ence between the flavour states equals n/2. This can be considered as the definition of 
refraction width, which by eq. (0) can be written in the general form: 

do = ^, (14) 

where V m is the neutrino-medium potential and n m is the concentration of the relevant 
scatterers in the medium. 

The minimum width, d m i n , is inversely proportional to tan 26, which represents prop- 
erties (the mixing) of the neutrino system itself. The smaller the mixing 9, the larger is 
the width d min needed for strong transition. 



2 It can be checked that different choices of the condition M) lead to analogous results. For instance, 
taking P Ue ^, Vft > | we find 

4 fi = -do = = 5.41 • 10 8 A cm- 2 , 

° 3 3V2Gf 



and, for P Ve ->v u > 



4 /4 = -d = = 2.7 • 10 8 Acm- 2 . 

3 3V2G F 
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The condition (f|) can be generalized. It corresponds to the case of initial state co- 
inciding with a pure flavour state. In general one can require that the change of the 
probability to detect a given flavour a is larger than 1/2: 

AP = P s {v a ) - Pi(v a ) > ~ , (15) 

where Pi and Pf are the initial and final probabilities. The condition ((^) corresponds to 
Pi{ v ij) = 0) so that Pf(v a ) = Py^u . Taking Pj = 1/4 and Pf = 3/4, we get in a similar 
way: 



3 mm 3 v / 2G F tan2fl 



d\/2 — -^rairi ~ ~ — ; ~ • (16) 



This d\/2 is the extreme value, however for most practical situations the condition (^) is 
more relevant, and from here on we will use the the width d m i n determined in @. 

In what follows we will show that for all the other density profiles the width din 
required by the condition (||) is larger than d min . 

2.2 Uniform medium with density out of resonance 

For n e ^ n r e es the inequality (|5|) can be satisfied only if sin 2 29 m > |, which means that the 
density is required to be in the resonance interval: n r e es (l — tan 20) < n e < n r e es (l + t&n26). 
At the edges of the interval we get the width 

di/a = (-hr Q ± V2d mn , (17) 



2G F Vtan20 

which is larger than d m i n . For other values of the density in the resonance interval we 
have d m i n 

2.3 Medium with varying density 

In general the neutrino propagation has a character of interplay of resonance conversion 
and oscillations. Two conditions are needed for strong transition: 

1) Resonance condition: the neutrinos should cross the layer with resonance density. 

2) Adiabaticity condition: the density should vary slowly enough. This condition can be 
written in terms of the adiabaticity parameter 7 at resonance: 



7 < 1 (1 
dn 



2E cos 29 

7 



Am 2 sin 28 n e 



dL 



(19) 



3 It can be checked that the width di/ 2 in eq. ( p7[ ) is larger than d m i n for small mixing: sin 29 <, 0.3. 
We will use this condition as criterion of smallncss of the mixing. 
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Notice that both the conditions 1) and 2) are local, and can be fulfilled for arbitrarily small 
widths of the medium. Clearly, they are not sufficient to assure a significant conversion, 
and a third condition of large enough matter width is needed. 

Let us consider a linear density profile with length 2L and average density equal to 



the resonance one, so that n r , 



+ An and n r 



n. 



An. Denoting 9\ m and 



#2m the mixings in the initial and final points, we find that in the first order of adiabatic 
perturbation theory the conversion probability is given by: 



where 



1 1 



cos 29 -i „, cos 29- 



2m 



— - sin 2Q\ m sin 2#2 m cos 



7 



-/(*) 



— 2 sin(2#i m — 292 m )a(x) cos 



(20) 



X 

/(*) 

a(x) 



27T7" 



ln(x + \/l + x 2 ) + xy/T 



x z 



dy 



o l + y 2 



cos 



20 



f(y) 



(21) 



For 7 < 1 we get from eqs. ([H]) and (|2lD : 



d\/2 — d-mir. 



1+1 



7T 



7 



(22) 



dmiry and for weak violation of 



This expression shows that for the adiabatic case di/2 
adiabaticity the minimum width increases quadratically with 7. We remark that in this 
case the effect is dominated by oscillations with large (close to maximal) depth. The 
change of density gives only small corrections. 

Let us consider now a situation in which the resonance adiabatic conversion is the main 
mechanism of flavour transition. A pure conversion effect is realized if the initial neutrino 
state that enters the medium coincides with one of the eigenstates of the Hamiltonian 
in matter, and the propagation in matter is adiabatic. In this case no phase effect, and 
therefore no oscillations occur. Let us denote rij and n/ the initial and final densities of the 
medium, and suppose the initial state is z/j = z/ 2m = sinOmVe + cos 6^1^. The probability 
to find a in this state is -Pj(^) = cos 2 6 m {ni). The state evolves following the change of 
density, so that it remains an eigenstate of the Hamiltonian, and the probability to find 
Up in the final state is P/(z/ M ) = cos 2 9 m (rif). Since the initial state u^ does not coincide 
with a pure flavour state we will use the condition ( |T5D as criterion of strong matter effect. 
Inserting P; and Pf in ([15]), we get the condition for di/ 2 : 



cos26 m {rif) 



cos 29 r 



[nil 



(23) 
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Taking the initial and final values of the density as n, = n r e es + An and n/ = n r e es — An 
(An > 0), and using the definition (f|) we find that the equality (^) leads to 

An = < es -J=tan2# . (24) 

Clearly, for a given An the size of the layer, and therefore its width, depend on the 
gradient of the density which can be expressed in terms of the adiabaticity parameter 7, 
eq. (0). We get: 

2Ecos26 1 , 

n e {L)dL = - — - , 2o/1 -an e , (25) 
Am 2 snr26>7 

and then integrating this equation we obtain: 

2Pcos2fl 1 

d=- — ^— -An . 26 

Am 2 sin 2 26 7 v ; 

Finally, inserting the expressions fl24|) and (|]) in eq. (|26"D we find: 

4 1 

di/2 — — j=—d m i n . (27) 



7T 



Let us comment on this result. As far as the adiabaticity condition is satisfied, the change 
of probability does not depend on the density distribution; it is a function of the initial 
and final densities only. If An is fixed, the decrease of the width means the decrease 
of the length L of the layer, and therefore increase of the gradient of the density. This 
will lead eventually to violation of the adiabaticity condition. Thus, the minimal width 
corresponds to the maximal 7 for which the adiabaticity is not broken substantially. 

For strong adiabaticity violation an increase of di/ 2 is expected, due to the increase 
of the minimum An required by the condition (|I5|) , and therefore of the correspond- 
ing length. This can be seen if we consider the previous argument taking into account 
the effect of the adiabaticity breaking from the beginning. Using the Landau-Zener level 
crossing probability Plz = exp(— 71/27), which describes the transition between two eigen- 
states, we get, instead of (p3|): 

(1 - 2P LZ )(cos20 m (n/) - cos2ft m (n i )) = 1 , (28) 

where we have averaged out the interference terms. Then instead of eq. (|24j) we get 



An = n r e es . 1 tan 26 . (29) 

'l6Pl z - 16P LZ + 3 



Finally, the condition for d 1 / 2 can be written as: 



4 

d±/2 = ; ^=d m i n . (30) 

Tt^JWPlz - 16P LZ + 3 
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For 7 — > eq. (^) gives <ii/ 2 — > oo, according to the fact that the density changes 
very slowly and therefore the width needed to have significant conversion increases. With 
the increase of 7 the width d\/2 decreases and has a minimum at 7 ~0.7, for which 
we find d\/2 — 1.5d min . With further increase of 7 (7 <: 0.7) the width d\/2 increases 
rapidly. According to (j30|) it diverges for P — > 1/4, when 7 — > 7r/(41n2) ~ 1.13. This 
value corresponds to the case in which the adiabaticity violation is so strong that even 
an infinite amount of matter is not enough to satisfy the condition (|T5|). Thus, we have 
found that also in this case g?i/ 2 > d m i n . 



2.4 Step-like profile 

As an extreme case of strong adiabaticity violation, let us consider the profile consisting 
of two layers of matter, having densities n\ = n r e es + An and n 2 = n r e es — An (An > 0), 
and equal lengths L\ = L 2 = L. At the border between the layers the density has a jump 
of size 2 An. We fix L = l res /8, so that d = d min . The result for the conversion probability 
can be computed exactly: 



pstep 



2 • 2 
S Sill 



7T 



2 2 

s c 



( 11 

1 — cos — 
V4s 



(31) 



where we denote the mixing parameters in the two layers as sin2# 2m = sin26 l lm = s, 
cos26 l 2 m = — cos26 l lm = c. In absence of the step (An = 0), P£^„ equals 1/2, recovering 
the case n e = n r e es = const. The probability ([H]) decreases monotonically as An increases. 
Expanding in 5 = (An/n r e es tan 28) 2 we get^: 



pstep 



1 



7T 



0.02 5 . 



(32) 



According to (|32|) , for d = d we have P*^* v < 1/2. This implies that, to have P^ v 
1/2 one needs d\/2 > d min . 



2.5 Castle-wall profile 

The profile consists of a periodical sequence of alternate layers of matter, having two 
different densities n\ and n 2 . We denote the corresponding mixing angles as Q\ m and 62m- 
In this case, a strong transformation requires certain conditions on the oscillation phases 
acquired by neutrinos in the layers [JIT}| ; therefore the transformation is a consequence 
of the specific density profile, rather than of an enhancement of the mixing. Suppose 
n% n r e es and n 2 = 0, and take the width of each layer to be equal to half oscillation 
length, so that the oscillation phase acquired in each layer is 7r. It can be shown that 
for small 9 this is the condition under which the conversion probability increases most 

4 This approximation proves to be very good (relative error < 0.5%) for < S < 1, i.e. for ni and n 2 
in the resonance interval. 
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rapidly with the distance. As a function of the number iV of periods (a period corresponds 
to two layers), the probability is given by 

P Ue ^(N) = sm 2 (2NA6) , (33) 

where AO = Q\ m — $2 m = Q\ m — &■ Using the approximation 2$i m — 2$ ~ sin 28i m — sin 28, 
and expanding sin26*i m in rii, we get: 

" 1/2 = " ndmm ■ (34) 

Again, we find that di/ 2 > d m i n . 

Thus, for all the known mechanisms of matter enhancement of flavour transition (res- 
onant oscillations, adiabatic conversion, parametric effects), we have found that the width 
d\j2 is larger than d m i n , which is realized for the case of uniform medium with resonance 
density. In fact the constant profile with resonance density could be expected from the 
beginning to represent an extreme case: this profile is singled out, since it is the simplest 
distribution with the density fixed at the unique value n r e es . 

It is worthwile to introduce also the total nucleon width. Let us consider a medium 
made of electrons, protons and neutrons with number densities n e , n p and n n . Defining 
the number of electrons per nucleon as Y e = n e /(n n + n p ), we can write the total nucleon 
width that corresponds to d Q as: 

doN = 77 • (35) 

* e 

We can also introduce the total mass width d p : 

dp = m N d 0N = m ^ d ° . (36) 

For electrically and isotopically neutral medium (n e = n n = n p ) , eq. (|35|) gives: 

d 0N = 2d , (37) 

and numerically: 

d 0N = 4.9 • 10 32 cm" 2 , d p = 8.16-10 8 g-cm- 2 . (38) 

3 Generalizations 

In this section we generalize the previous results to active-sterile transition and to the 
case of flavour- changing induced conversion. We also discuss the small width limits. 
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3.1 Active-sterile conversion 



In this case the scattering both on electrons and on nucleons contributes to the conversion, 
and the effective potential for an electron neutrino in an electrically neutral medium equals 



V = V2G F n e (l 



n, 



2n r 



(39) 



Thus, the results for v e — v s transition can be obtained from those for v e — by the 
replacement n e — ► n e (l — n n /2n p ). For the refraction width we get immediately: 



d (u e 



dn 



2n r 



2flp Ti n 



In particular, for an isotopically neutral medium eq. 

d (u e -> u s ) = 2d . 



gives 



(40) 



(41) 



Notice that, for highly neutronized media (n n ^> n p ) the width d (v e — > v s ) gets signif- 
icantly smaller than do- In this case, however, the physical situation is more properly 
described by the total nucleon width d N defined in eq. (j35|), since the effect is mainly 
due to the scattering on neutrons. We find: 



which gives in the limit n n —>■ oo: 



2d 



Tip "T Tl n 



2flp TL n 



doN^e 



2d, 



o • 



(42) 



(43) 



similarly to eq. (41). 

For the — v s case, the potential, and consequently the width, can be obtained by 
replacing n e — > n e (—n n /2n p ), which gives 



do(v„ 



2n p 
do- 



11,. 



(44) 



For isotopically neutral medium, eq. fl44|) reduces to eq. (|41|). For highly neutronized 
media, the argument is analogous to the one for the v e — v s case, and we get the same 
result as in eq. (f43"D. 



3.2 Oscillations induced by flavour changing (FC) neutrino-matter 
interactions 

In this case the neutrino masses can be zero, or negligible, and the flavour transition is a 
pure matter effect. The Hamiltonian of the system has the following form M: 



H = V2G F (° 

\en f en f J 



(45) 
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where rif is the effective number density of the scatterers, and e and e' are parameters of 
the interaction. 

As follows from eq. (|45|), in a uniform medium the neutrinos oscillate with transition 
probability: 

^W = i ^sin 2 ( 7 r y ), (46) 
/ = _JE^ L_ (47) 

We assume e' < e, which is needed to have a significant oscillation amplitude. Using eq. 
(D, we get: 

d FC * I He 

1/2 2y/2G F V4e 2 + e /2 n f 

dn n e . 
= / — • 48 



Notice that the factor (n e /n/) tan2#/\/4e 2 + e' 2 implies that d\9 % can be significantly 
smaller than d m i n , and the oscillation effect can be observed in media of smaller width. 
For a FC neutrino interaction with up (or down) quarks and isotopically neutral medium 
we have n e /rif = 1/3, and therefore: 

d\% ~ ld • (49) 
' be 

Notice that there are two sources of decrease of the width: the factor 2 is given by the 
presence of two off diagonal terms in the Hamiltonian fl45f) and the factor 3 is due to the 
larger number of scatterers. Taking e = 1, we find the value: 

d{f 2 ~ 1.36- 10 8 g-cm~ 2 . (50) 

For a density n = 4g • cm -3 (Earth's crust), this corresponds to the distance L = 337 Km, 
which is comparable to the length of the present long base-line neutrino experiments: K2K 
(base-line 250 Km) and ICANOE (740 Km). 

3.3 The small width limits 

In a number of situations (see section 5) the width of the medium is smaller, or much 
smaller, than d m i n . We consider, then, the matter effect on oscillations in the limit 
d/d m in 1. Introducing the two variables: 

A = — = Ue (51) 

- [res/4 P ~ n r e es ' 
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we can write the small width condition as: 

d/dmin = Ap < 1. (52) 
We focus on two important realizations of this inequality: 

1) Small size of the layer and density close to resonance. As we have shown in section 2, 
a strong transition requires n e ~ n r e es . In case of small width this implies small size of 
the layer. Therefore we have A 1 and p ~ 1. In this case d/d m i n ~ A. We expand the 
oscillation probability (^]) in A at the lowest (nonzero) order, and find that the matter 
effect vanishes quadratically with d/d min : 




2) Small density of the medium and length close to the minimum value / res /4. Another 
condition of strong conversion is to have the size of the layer of the order of the oscillation 
length. According to eq. (|52"D, this means that the density is small. Thus, we have p <C 1 
and A ~ 1, and therefore d/d m i n ~ p. In order to give a phase-independent description of 
the matter effect, we perform an expansion in p of the oscillation amplitude: 

sin 2 26 m - sin 2 26 = 2p sin 2 26 cos 26 ~ 2-— sin 2 26 cos 26 . (54) 

Unlike the previous case, the relative matter effect is linear in d/d m i n . 



4 Refraction of high energy neutrinos 

In this section we examine the refraction of high energy neutrinos (s ^ M^), both in 
matter and in neutrino background. 



4.1 High energy neutrinos in matter 

Let us consider the propagation of high energy neutrinos in medium composed of protons, 
neutrons and electrons. The expressions (H) and (|II]) refer to the low energy range, 
s <C M^y, where s is the center of mass energy squared of the incoming neutrino and the 
target electron, and My/ is the mass of the W boson. The general formulas, valid for high 
energies too, can be obtained by restoring the effect of the complete propagator of the W 
boson in the expression of the potential (an analogous argument holds for the Z boson). 

Let us consider v e — conversion. Since the refraction effects are determined by the 
real part of the propagator, the potential ([1]) is generalized as: 

V = V2G F nJ(q 2 w ) (55) 
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where q w = q 2 /M w and = Fw/Mw] q and IV are the four momentum and the width 
of the W boson. 

The only contribution to the potential ( |55D is given by the forward charged current 
scattering on electrons (w-channel exchange of W), for which q 2 ~ — s. Therefore, in- 



troducing sw = s/M$y, we have q w ~ —sw 
refraction width do using the definition (PI), 
find: 

a i \ ! ^° 



From the potential ( pq) we can find the 
For the nucleon refraction width 



do 

Ye 



1 + s w ) 



we 



(57) 



where we have neglected the width jw Eq. fl5"7|) shows that g?o7v( s w0, and therefore 
^mm(siy), increase linearly with sw above the threshold of the W boson production. 

For the active-sterile conversion, one has to take into account also the neutral current 
interaction channel (t-channel exchange of Z), for which q 2 = 0, so that the low energy 
formulas ([H-jlOD are still valid. For — v s only neutral current processes are involved, thus 
the low energy result, eq. ([44D , holds at high energies too. In contrast, for the v e — v s case 
both charged and neutral current interactions contribute, and for an electrically neutral 
medium the high energy potential can be written as: 



V = V2G F n e 



71, , n 



1 + sw 2rir 



(58) 



The second term in eq. flog ) does not depend on sw, thus coinciding with the correspond- 
ing term in the low energy expression (|39|). The potential fl58|) gives the nucleon refraction 
width: 

1 + s w 



doN{Sw) — 2c?o 
For isotopically neutral medium (Y e = 



(3Y e - 1) - s w ( 
1/2) we get: 



y: 



d N{sw) = 4c? 



1 + Sw 



1 - s w 



(59) 



(60) 



The width d N(sw) diverges for sw —> 1 (see fig. |I|). 

At high energies inelastic interactions and absorption become important: at sw ~ 1 
the imaginary part of the interaction amplitude is comparable with the real part. In fig. 
HI we show the refraction width d p = m^doN for u e — and u e — v s conversion and the 
absorption width d^ |15, 16| as functions of the neutrino energy E in the rest frame of 
the matter. We have considered isotopically neutral medium, Y e = 1/2. The absorption 
width d a b s is dominated by the contribution of neutrino-nucleon scattering; it decreases 
monotonically with the energy E. In contrast, d p starts to increase at sw ~ 1, which 



corresponds to E = 10 6 -=-10 7 GeV, according to eq. (|57D . For E ~ 10 6 GeV absorption and 
refraction become comparable; at higher energies, the former effect dominates: d a b s ^ d - 
This means that for a v e with energy E > 10 6 GeV the conversion in matter is damped 
by inelastic interactions and absorption JT7j^H|, therefore a strong conversion effect can 
not be expected. 
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Notice that for small mixing angle 6 the minimum width d m i n is significantly larger 
than the refraction width do, therefore the absorption starts to be important at lower 
energies. Taking, for instance, sin 26 = 0.3 we have d m i n ~ do/ sin 26 ~ 3.3do, and find 
that d a b s & d min already for E ^ 5 • 10 5 GeV. 




Figure 1: The width d p = m^doN for v e — (dashed line) and for v e — v s (dotted line) 
channels, and the absorption width, c? o6s , for the electron neutrino (solid line), as functions 
of the neutrino energy. We have considered isotopically neutral medium, Y e = 0.5. The 
data for d a b s are taken from ref . [|T3| . 

Let us consider now the matter effect for conversion of antineutrinos. For u e — 
channel the only contributing interaction is the u e — e scattering with W exchanged in 
the s-channel. In this case q 2 = s, and using eq. (|5lf) we get: 

dmM = ■ (61) 

This function has a pole at sw = 1, i-e., in the V^-boson resonance. The pole appears 
because the amplitude becomes purely imaginary in the resonance, so that the potential 
is zero. The width doN^sw) diverges for sw —> oo, due to the l/sw decrease of the 
amplitude. The function (|6T|) has two minima: 

do N {s™ n ) = 2 lw ^ = 2 lw d 0N {s w = 0) at s$ n = l± lw . (62) 

I P. 
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Numerically doN(s\v n ) = 0.05 doN(sw = 0), which shows that refraction effects are en- 
hanced close to the W resonance. However, in this region inelastic interactions become 
already important. 

For v e — v s channel the contribution of neutrino-nucleon scattering should be included, 
and for electrically neutral medium we find: 



^07v( s iy) — 2g?o 



;i-sw/) 2 + 7w/ 



(3Y e - 1) + 2s w (l - 2Y e ) - s 2 w {\ - Y e ) - 7^(1 - Y e ) 
In the case of isotopically neutral matter eq. fl63f) gives: 

;i -swf+iw 



(63) 



dw(sw) = 4c? 



1 Syy 7jy 



(64) 



which has the value Ado in the limits sw "C 1 and sw 3> 1, and a pole at sw — 1- Similarly 
to eq. ( j62| ) we find the minima: 



^07v(sw in ) = 47^/^0 = 7wd 0N (s w = 0) at = l±j w . 



(65) 



In fig. we show the refraction width oL = mjsidoN for z/ e 



• and z/ e • 



v s channels and 

the absorption width for the electron antineutrino, d a b s ||16|| , as functions of the neutrino 
energy. We have considered isotopically neutral medium. For energies outside the W 
boson resonance interval the main contribution to d a b s is given by the neutrino-nucleon 
scattering; at Sw — 1 the effect of the resonant u e — e scattering dominates, providing 
the narrow peak. It appears that absorption prevails over refraction {d a bs < do) for 
E ^ 6 • 10 6 GeV, corresponding to d p ~ 6 ■ 10 7 g ■ cm -2 , for both v e — and u e — v s cases. 

The effect of absorption on neutrino conversion starts to be important at lower ener- 
gies: for sin 29 = 0.3 we find that d a b s ^ d min at E ^ 6 ■ 10 5 GeV. 



4.2 High energy neutrinos in neutrino environment 



Let us consider a beam of neutrinos which propagates in a background made of neutrinos 
of very low energies^. This could be the case of beams of low energy neutrinos from 
supernovae, or high energy neutrinos from AGN and GRBs, or neutrinos produced by 
the annihilation of superheavy relics, etc.. We assume that the background consists of 
neutrinos and antineutrinos of various flavours, with number densities rij (i = u e , u e , u^, 
etc.). In the case of relativistic neutrino background we assume its isotropy. 

The potential for a neutrino u a (a = e, fi, r) due to neutrino-neutrino interaction can 
be written as: 



V Ua (s z ) = V2G F 



n 



>J{sz) -n Da f(s z ) + ( n i~ 



m 



t = V e M a M T 



(66) 



5 We will not consider the conversion of neutrinos in the background itself, which can significantly 
affect the flavour content of the background. 
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Figure 2: The width d p = mudoN for v e — (dashed line) and for v e — v s (dotted 
line) channels, and the absorption width, d a b s , for the electron antineutrino (solid line), 
as functions of the neutrino energy. We have considered isotopically neutral medium, 
Y P = 0.5. The data for gL/,, are taken from ref. |T6l . 



where the propagator function f(sz) has been defined in eq. (f)6|). Here sz = s/M| and 
'-/z = Tz/Mz; Mz and Tz are the mass and width of the Z-boson. The first term in eq. 
(|66|) is due to v a — v a scattering with Z-boson exchange in w-channel, and the second 
term is the contribution from v a — u a annihilation. 



4.3 Neutrino conversion in CP-asymmetric neutrino background 

As a first case we consider a strongly CP-asymmetric neutrino background, and suppose 
Hi ^> rij, so that we can neglect the contributions of antineutrinos in (0). For simplicity, 
we assume equal concentrations for the neutrino species: n Ue = n u = n Vr . In terms of 
the total number density of neutrinos, n u = n Ue +n u +u Ut , the potential ( j66|) reduces to: 



V Ua (s z ) = V2G F n u 
The potential for the antineutrino is given by V Pa (sz) = —V Ua (—s z ). 



(67) 
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Let us now find the refraction width do and d m i n for various channels. 

1). For the active-sterile conversion, v a — i/ s , the potential (^) coincides with the 
difference of the potentials for the two species, and therefore, by eq. flT4|) , it gives imme- 
diately the refraction width of neutrinos: 



d (sz) = d 



l + lf(sz) 



The width d Q (sz) is constant for s z C 1 and sz 1: d (sz <C 1' 
10 32 cm~ 2 , and d (s z » 1) ~ c£ (see fig.||). 



(68) 

3d /4 = 1.84 • 



Let us now compare the refraction and absorption effects. The main contribution 

| is given by the v a — v a and v a — vp {(3 ^ a) scat- 
decreases monotonically with sz and at sz 1 it takes the 



to the absorption width d a b s 
terings. The width d abs 



value d abs (sz 3> 1) — n / (2G 2 F M^) ~ 3.6 • 10 33 cm -2 . Due to its non-resonant be- 
haviour, d a b s is larger than do for any energy of the neutrinos: at Sz 3> 1 we find that 
do/ dabs — GfM\/\[2 = 7raw/(2 cos 2 9w) — 0.1, where 6w and aw = 5 ,2 /47r are the weak 
mixing angle and coupling constant. Therefore, d a b s is also larger than the minimum 
width, d min , for sin 29 k d /d a b s — 0.1. 



2). For the conversion of an active antineutrino into a sterile species, v a — u s , we get 
the width: 

1 



do(sz) = d, 



' + -f{sz) 



(69) 



This function (see fig.|3|) has a resonant behaviour with minima at sz — l±7z: do{l— jz) — 
(1/67^ + l)" 1 ^ — d /7 and d (l + j z ) — (I/672 — l) _1 rf — <i /5. Outside the Z-boson 
resonance d (sz) is constant: d (sz < 1) - 3d /4 and d (sz ^> 1) — do- In the range 
sz ~ 1 inelastic scattering and absorption become important. We evaluate the absorption 
width dabs for antineutrino in neutrino background using the plots in ref. p0| . For sz < 1, 
the width d a b s decreases with the increasing sz', at sz — 1 it shows the characteristic 
peak due to the resonant u a — v a scattering. For sz k 1 the absorption width increases 
with sz up to the limit d a b s { 
scatterings (/? 7^ a) and v a - 
for s z ^ 0.8 
cWsz = 0.8) 



s z ^ 1) ~ 3 • 10 cm , due to the contributions of v$ 



^ 7 



interaction in the t-channel. We find that d a b s ~ do(sz) 
10 3 GeV 2 ) and for s z k 2 (s k 1.7 ■ 10 4 GeV 2 ). Furthermore, 



10 32 cm~ 2 and d (sz 



2) ~ 4 ■ 10 32 cm- 
get that d a b s ~ d m { n for s z 0.7 and s z 

^ 2.4, drain 



Taking sin 29 = 0.3 we 
0.7) ~ 3 ■ 10 32 cm" 2 and 
dmin(s z = 2.4) ~ 10 33 cm- 2 . 

Notice that, in contrast with the conversion in matter (see figs.|l] and |2|), for neutrinos 
and antineutrinos in neutrino environment we can have d a b s k d m i n (sz) even in the high 
energy range, Sz k 1: in particular, we find that d a b s (sz 3> 1) *J d m i n (sz 3> 1) for 
sin 29 k 0.1. 



3). Let us now consider the active- active conversion, v a — vp. Assuming equal concen- 

u v = n UT , we find from eq. (|67j ) that the 



trations of neutrinos of different flavours, n v 
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Figure 3: The dependence of the refraction width d for v a — v s (dashed line) and u a — u s 
(solid line) channels on sz in strongly CP-asymmetric background, n$ ^> rq. Equal 
concentrations are assumed for the various flavours. 



difference of the potentials of the two species equals: 



V Va - V u = -V2G F n v \f{ 



/(- 



(70) 



where s l z = s*/M| (i = a, /?), and s l is the center of mass energy squared of the incoming 
and the background neutrino of the same type i. If the species v a and up in the background 
have different energies we get that s% ^ s§, and therefore AV a ^ ^ 0, leading to matter 
induced neutrino conversion even if v a and vp have equal concentrations. This situation is 
realized if the background neutrinos v a and vp have different masses, e.g. m Ua > m U g, and 
are non-relativistic. Denoting by E the energy of the neutrino beam, in the rest frame 
of the background we have s l = 2miE, and thus s^/s§ = m Uct /m U0 > 1. The condition 
s z 7^ s z i s achieved also if one of the neutrino species is relativistic and the other is 
not: m Va ^> Ep^> m v „, where Ep is the energy of vp in the background. Assuming the 
isotropy of the neutrino gas, we have that ~ 2EpE. 
Using (|70|) and (0), we get the refraction width: 



do(4) = 3^o /(-4) " f(~4) 



3d 



;i + S f)(i + 4) 



l S Z S Z/ 



(71) 
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The function (|n]) diverges for s l z 
limit, sz <C 1, it reduces to do(s z <C 1) 
the realistic case s z ^ 1 and s z <C 1, eq. 

do(al) = 



oo and — > 0. In particular, in the low energy 



3dn 



3do/(2EAm), where Am 
H) can be written as: 



m 



For 



(72) 



which approaches the minimum value 3g?o when s z 3> 1 (see fig.[|). Taking the maximal 
realistic values for the mass and energy of the neutrino, m Ua = 5 eV and E = 10 22 eV we 
get s z ~ 12 at most, so that do(s z ) ~ 3.5do- 



4). For the v a — vp channel the effective potential equals: 



and therefore we get the width: 

d (s z ) = 3d 



f{s a z) ~ /(4) 



(73) 



(74) 



Due to the resonant character of the function f(s z ), the width d (s z ) has the following 
features (see fig||): 

(i) It reaches the local minimum do(s z ) ~ 67^ do ~ do/6 when one of the s z 's is at 
resonance and the other is outside the resonance: s z ~ 1 and s§ 7^ 1 (or vice versa). 

(ii) The absolute minimum d (s l z ) ~ 37^d is achieved when s z ~ 1 + 7z and s§ ~ 1 — 7^ 
(or vice versa). These conditions can be satisfied for certain relations between the masses 
of the background neutrinos. For non-relativistic background: m Va /m Uf3 = (1 + 7z)/(l — 
Iz) ■ 

Notice that for s l z discussed in (i) and (ii) the effects of inelastic scattering and ab- 
sorption can be important. 

(iii) If one of the s l z s is far below the resonance and the other is far above (e.g. s z ^> 1 
and s z <C 1) then do(s z ) ~ 3do. 

(iv) d (s z ) ^> d if both the s % z s are far below or far above the resonance. 
Obviously, for strong CP-asymmetric background with ^> rij the results for v and 

v channels should be interchanged. 



4.4 Neutrino conversion in CP-symmetric neutrino background 

Let us now consider the neutrino conversion in a CP-symmetric neutrino background, 
rii = nj, with n Ue = n u = n Vr . In this case the potential (|66|) can be written as: 

V Va = V2G F n Va if(-s z ) - f{s z )} . (75) 
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2.5 



Figure 4: The dependence of the refraction width do for v a — up (dashed line) and V a — up 
(solid line) channels on sz in strongly CP-asymmetric background, rij ^> nj. Equal 
concentrations are assumed for the various flavours. We have considered ^> s§ ~ 0. 



It vanishes in the low energy limit sz — > 0, but it is unsuppressed at high energies, sz *S 1, 
leading to significant matter effect. We will consider the conversion of neutrinos; due to 
the CP-symmetry antineutrinos undergo analogous effects. 

1). For u a — v s conversion from the potential (|75|) we find the refraction width for the 
neutrinos of flavour a: 



do(sz) 



do\f( 
d 



sz 



f{sz) 



i + 4 + 7l) 2 -44 



2s z [(l 



+ 7l] 



(76) 



For sz ^ 1 it behaves as do(sz) — do/{2sz) and for > 1 we have d Q (sz) — d sz/2 (see 
fig.^|). The width (|76|) has two minima: 



d, 



2-fzdo 



at S' 



1 ±7 Z . 



(77) 



Numerically, d (sf in ) = 0.055 d ~ 1.35 • 10 31 cm~ 2 . 

The absorption width d a b s is dominated by u a — v a annihilation, with a resonance peak 
at sz ~ 1. Using the results of ref. |2D| we find that d a b s is larger than d (sz) outside 
the Z-boson resonance, and the two quantities are comparable at sz ~ 1 or at sz ^> 1. 
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For sin 26 = 0.3 we find that the minimum width d m i n is larger than d a \, s in the range 
0.7 £ s z £ 1.6, corresponding to 6 • 10 3 GeV 2 £ s £ 1.3 • 10 4 GeV 2 . At the edges of this 
interval the width d min takes the value d min (sz = 0.7) ~ d min (s z = 1.6) ~ 3 • 10 32 cm~ 2 . 

2). For the v a — v@ channel, eq. ( [75|) gives the difference of potentials: 

= V2G F n Va {[f{s a z ) - m)) - [f (-4) - /(*£)]} • (78) 
The corresponding refraction width equals: 



do{s % z ) = d 



[/(-«!) -/(«!)] -[/(-4)-/(4)] * • (79) 



We find that c?o( s z) ^ <^o when both s% and s§ are outside the Z-boson resonance and 
do(s l z ) takes its minimum values when either s% or s z is close to the Z-resonance: 

If s z ~ 1 + '-fz and s z ~ 1 — 7^ (or vice versa) g?o(s^) has the absolute minimum 
e?o(s^) — 7z^o- F° r s z ~ 1 an d s z 7^ 1 ( or y i ce versa ) we have the local minimum 
do(s l z ) ~ 2 , -fzdo- Notice that in the realistic case s% ^> s z ~ the width ( |79|) reduces 
essentially to the one in eq. (0). At resonance, where do(s l z ) has minima, the effects of 
inelastic collisions and absorption are important. 



5 Applications 

The results derived in the previous sections are now applied to some physical situations 
of interest. After a brief discussion of well known cases, like the Earth, the Sun and 
supernovae, we present results for neutrinos in some new astrophysical environments. We 
find that significant matter induced conversion can be expected for neutrinos crossing the 
dark matter halos of clusters of galaxies and for neutrinos from cosmologically distant 
sources. 



5.1 Minimum width condition and bounds on the mixing 

As follows from the analysis in section 2, a significant neutrino conversion in matter 
requires the fulfilment of the minimum width condition^: 

d ~ dmm = SlO ■ (80) 

This condition is independent of the density distribution, and therefore of the specific 
matter effect involved. Thus the knowledge of the width d allows one to conclude about 
the significance of the matter effect even if the density profile is unknown. This is the 



6 This condition refers to the requirement of conversion probability larger than 1/2, eq. (||). In some 
circumstances, however, even a small effect, with conversion probability Pcl/2 can be important. 
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Figure 5: The dependence of the refraction width do for v a — u s and v a — v s channels on 
Sz in CP-symmetric neutrino background. 

case of some astrophysical objects for which estimates or bounds on d can be obtained 
directly by observational data with no assumption on their internal structure. 

In the Table [l] we show the parameters of interest of some objects, toghether with the 
values of the ratio 

d 

For r < 1, and small mixing angle, the condition ( |5DD can not be satisfied, thus no 
significant neutrino conversion is expected. Conversely, for r ^ 1, (|80|) can be fulfilled 
and gives the bound on the mixing: 

sm26Z- = ^. (82) 
r d 

Notice that our analysis holds for small mixings: sin26 l <C 1. For applications we assume 
sin 28 ^ 0.3, for which we find from eq. (jS2] ) that the minimum width condition is satisfied 
for r ^ 3. 

The inequality (|82|) can be considered as the sensitivity limit for the mixing angle that 
can be achieved by studies of neutrino conversion in a layer of given width d. The real 
sensitivity can be however much lower than the absolute limit given by the condition fl8~2"D . 
This is related to the fact that in the case of varying density only part of the total amount 
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of matter effectively contributes to the conversion. Introducing the corresponding width 
d conv we have the condition: 

sin2#>-^, (83) 

instead of the ( |§2] ) . 

Let us find the expression of d conv for a medium with monotonically varying density. 
As discussed in section 2.3, the transition occurs mainly in the resonance layer. Using the 
result (|24]) we get: 

d conv = n res 2An = n ™>i t an 29 , (84) 
e dn 

where l n = |(^)~ 1 |res ri e es - Inserting the expression fl84|) in the condition (p3|), we find: 

sin 2 29 > -^L . (85) 

Clearly, d conv could be much smaller than the total width d of the object, so that the 
condition (^5|) on the mixing could be much stronger than ([82]). Notice that the bound 
fl35p is quadratic in sin 29. Using the definition (|l~9|) of the adiabaticity parameter, 7, 
the condition ( p5|) can be written as 7 < 4/(7rv^3), which corresponds to the adiabaticity 
condition close to its limit of validity. 

Another important issue is that the maximal sensitivity for the mixing 9 can be 
achieved for particular values of Am 2 /E, which depend on the specific density profile. 
As follows from (|8lS|), for constant (or slowly varying with the distance) l n the smallest 
sin 2 29 corresponds to the largest n™ s , and therefore to the largest values of Am 2 /E. 
This is the case of exponential density profile. For power-law profile, n e ~ L~ k , we get 
\l n \ = L/k, so that sin 2 29 ~ L . Taking k > 1, fulfilled by practically all the realistic 
profiles, we find that the smallest 9 is achieved for the smallest L, and consequently the 
highest values of n e and Am 2 /E. 

Notice that d conv is a local property which depends on the derivative in l n . Of course, 
the description given by d conv is not correct when the density profile is close to the constant 
one, so that l n —>■ 00. In this case d conv can be even larger than the total width d. Thus, 
the correct condition on the mixing can be written as: 

sin 20 > , r ^° ; r. (86) 

mm [d, d conv ] 



5.2 The Sun, the Earth, the Moon and supernovae 

For neutrinos crossing the Earth we consider two types of trajectories, corresponding to 
different values of the zenith angle 9 Z . For cos 9 Z =1 neutrinos travel along the diameter 
of the Earth, crossing the core and the two layers of the mantle. We get r =13.6, and 
therefore according to ( |52"D we could expect significant matter conversion for sin 2 29 ^ 
5 • 1CT 3 . However this maximal sensitivity, which would be achieved for uniform density 
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object 


density (cm 3 ) 


size (cm) 


r = d/do 


Earth: 








cos 9 Z = 1 


2.6 • 10 24 


1.26 • 10 9 


13.6 


COS 0* = 0.81 


1.5 ■ 10 24 


10 9 


6.4 


Sun 


~ 7 • 10 24 


6.96 • 10 10 


2600 


Moon 


~ 10 24 


3.48 • 10 8 


1.4 


Supernova 


3 • 10 33 


10 v 


10 9 


Universe 








(n v = n D ) 


1.5 • 10 4 


10 27 


3-10- 2 


Universe 








{n v » Tip) 


~ 10 5 


10 27 


0.3 


Galactic halo 


^ 2 • 10 6 


3 • 10 23 


5-10- 2 


Cluster halo 


£ 5 • 10 7 


3 • 10 24 


10 


AGN 




d ~ 10 22 -T- 10 23 cm~ 2 


10~ 10 -T- 10- 9 


GRB 


10 10 -T- 10 12 


< 5 ■ 10 15 


< 10~ 5 



Table 1: The density, the size and the matter width in units of refraction width, r = d/do, 
for various physical objects. The values given for the densities are averaged along the 
trajectories of the neutrinos. We quote the number density of electrons for objects made 
of usual matter, and the concentration of the neutrino background for the halos and the 
universe. For the Earth the results are given for two trajectories with different zenith 
angle 9 Z . The results for the universe correspond to redshift z = 5 for the cases of v a — u s 
and v a — v s in CP-symmetric and strongly CP-asymmetric neutrino background with 
rj v ^ 1. 



distribution, is not realized for the Earth profile. For small mixing, the difference between 
the densities in the core and in the mantle is larger than the resonance interval. As a 
result, the oscillations are resonantly enhanced either in the mantle or in the core, and 
only one of the two parts effectively contributes to the effect. At the same time, for 
certain ranges of Am 2 /E, different from both the resonance values in the core and in the 
mantle, parametric enhancement of oscillations occurs. Numerical calculations || give 
sin 2 29 <: 2 ■ 10~ 2 as best sensitivity. 

For cos # z =0.81 the trajectory is tangential to the core, and therefore it represents the 
path of maximal length in the mantle. In this case we find r ~6.4 and the sensitivity 
limit sin 2 29 2.5 ■ 10 -2 . Since this case realizes approximatively the optimal condition 
of uniform medium, we have good agreement with the results of exact calculations. 

In the case of the Moon, r =1.4, and therefore a large mixing is required: sin 2 29 ^0.5. 

A numerical integration of the density profile of the Sun [21] gives d ~ 1.5-10 12 g • cm -2 . 
Dividing this result by d p = m^doN, with Y e = 0.7, we find r ~2600. From the condition 
([82] ) we get then sin 2 29 ^ 1.5 ■ 10 -7 . This bound is remarkably weaker than the one 
obtained from the condition (jS5|): taking n r e es ~ 50v4cm -3 and l n ~ O.3_R , we get 
sin 2 29 ~ 2.4 • 10 -4 , in good agreement with the results of exact computations [@,|3|]. 

For supernovae the total width of the matter above the neutrinosphere gives r ~ 10 9 , 
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for which the condition (82) would lead to sin 2 29 ^10 18 . Using the density profile n e = 
n° e (R /R) 3 @, with R = 10 7 cm and n° e ~ 10 34 cm- 3 , from (§5|) we find sin 2 9 Z 10~ 8 , 
which agrees well with the results of numerical calculations ||. 

As shown in the previous examples, the maximal sensitivity for sin 2 29, given by the 
total width d, can be achieved in the case of uniform medium at Am? /E corresponding 
to the resonance density. Such a situation is realized for neutrinos crossing the mantle of 
the Earth. In the case of substantial deviations from the constant density, like in the Sun 
or in supernovae, the sensitivity is much lower. The stronger the deviation from constant 
density, the smaller d conv , and therefore the lower is the sensitivity. 

5.3 AGN and GRBs 

Let us now turn to high energy neutrinos from Active Galactic Nuclei (AGN) and Gamma 
Ray Bursters (GRBs) ||-g6j. 



In AGN, neutrinos are considered to be produced by the interaction of accelerated 
protons with a photon or proton background [f27) -|2~9|j . There is a hope that neutrinos from 
AGN with energies ^ 10 6 GeV could be detected by large scale underwater (ice) and EAS 
detectors [|30|j . 

The width of matter crossed by neutrinos in an AGN can be estimated on the basis 
of the existing data on the X-ray emission of these objects. The variability of the spectra 
suggests that the X-radiation is emitted very close to the AGN core ||24|| . The proton 
acceleration and therefore the neutrino production are supposed to happen in the same 
region. For this reason the width of matter crossed by neutrinos equals approximatively 
the one crossed by the X-radiation. For the later the experimental data [[H|-|33| give the 
value dAGN — (10~ 2 4 lCT 1 ) A cm -2 , therefore significant neutrino conversion in AGN is 
excludedF] (for a short discussion, see also ref. P^]). 

A rather successful description of the origin of GRB is provided by the fireball model 
|35l| , in which neutrino production is predicted to happen in an analogous way as in 



AGN [|T^,|36|]. A fireball can emit protons, detected as high-energy cosmic rays on Earth, 
accompanied by a flux of neutrinos. The requirement that the fireball should be trans- 
parent to protons gives an estimate of the width of the object: dcRB < d a t, s , where 
d a bs = 10-r- 100 A cm -2 is the total absorption width for the protons. It is possible to 
evaluate the width in a different way. An estimate of the electron number density in the 
fireball is given in ref. fL3 |: ugrb — (10 10 4 10 12 ) cm" 3 . Using this value, and taking the 
fireball mass in the range of star-like objects, M = (1 4 10) M Q , we can get the radius of 
the object, R GRB = 5 • (10 14 4 10 15 ) cm, and then the width: d GRB = 10 4 10 4 A cm" 2 . 
In agreement with the first argument, we see, then, that also in GRBs the matter effect 
on neutrino conversion is negligible. 



7 In the present discussion we have considered radial propagation of neutrinos from the inner to the 
external regions of the object. We have not considered neutrinos travelling through the core of the AGN. 
In this case a significant matter-induced conversion could occur, however neutrinos crossing the core are 
supposed to be a small fraction of the total neutrino flux produced. 
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5.4 Dark matter halos 



According to models |37| , part of the dark matter in halos of galaxies and clusters of 
galaxies should consist of neutrinos^]. Therefore neutrinos of extragalactic origin crossing 
the halo on the way to the Earth undergo refraction on the neutrino background. It 



was suggested in ref. |TJ] that, due to non uniform neutrino density distribution in the 
galactic halo, ultrahigh energy neutrinos can be resonantly converted into active and 
sterile species. 
Following ref. 



14 1 we consider a galactic halo composed of non-relativistic neutrinos 



and antineutrinos of the two species and v T . The electron neutrino is assumed to be 
lighter, and therefore less clustered, than and v T : n Ue <C rii, i = v^v T . We assume 
CP-symmetry of the background: n. 



nj. We take the density profile fH 
1 



n v vr) 



n 



o 



'l + (r/a)' 



(87) 



where n v = n u ^ + n„ T and a is the core radius of the halo. For galactic halos this radius 
is estimated to be a ~ 10 4- 100 kpc. An upper bound for n° is given by the Tremaine- 
Gunn condition P8fl : for identical fermions of mass m, the maximum number density 



allowed by the Pauli principle equals 



n ' 



6n 2 



(mv 



where v esc is the escape velocity of the particle: 

M x 



/2GMy 
V R ) 



540 



R 



1O 12 M 



30 Kpc 



Km/s. 



(89) 



Here M and R are the total mass and radius of the galaxy. In what follows we will take 
R ~ 3a, corresponding to the radius at which the density reduces at one tenth of its core 
value, n u (3a) = n°/10. From (|88|) and (|39^) we get: 



n ' 



1.7- 10 b 



m,, 



M 



,5evy v l012M ©/ v 10K p c , 

The integration of the profile (^) gives the matter width: 

d 



cm 



(90) 



71 o 
-n v a 



(91) 



Inserting the expression (|90| ) for n Tl 
d<8 - 10 28 



in (pi]), we find the upper bound for the width: 



m,, 



M 



5eV/ \10 12 M F 



10 Kpc 



cm 



(92) 



8 In what follows, we will not consider the heavy particles present in the halos, because their number 
density is much smaller than the one of neutrinos, although they provide the largest part of the mass 
of the halo. Furthermore, the amplitude of the forward scattering of neutrinos on neutrinos and on the 
heavy particles of dark matter are comparable, or the former is even larger. 
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According to ( p2\) the largest values of d are achieved for objects with big mass M and 
small radius a; so that compact halos represent the most favourable case. 

Let us now check the minimum width condition for — v s and — v e conversion in 
galactic halos. We use the refraction width d ( s ™ n ) — 0.055<io — 1-35 • 10 31 cm~ 2 given 
in eq. ([H]). This value was obtained for — v s conversion in CP-symmetric neutrino 
background, and is the absolute minimum of do(sz) (eq. ([7SD), realized at the Z-boson 
resonance, sz ~ 1. Notice that, under the assumption n Ue <C n,, the result rf (s™ n ) holds 
also for Up — u e conversion: due to the absence of electron neutrinos in the background, 
n Ve ~ 0, the potential (|75|) for v e is negligible, and therefore the electron neutrino behaves 
as a sterile species. 

From eq. ( P2"D we see that, for typical values of M and a of a galaxy, like for instance 
the Milky Way (M ~ 1O 12 M and a ~ lOkpc), the minimum width condition is not 
satisfied: d/d (sV} in ) ^ 5 • KT 3 . For the galaxy M87 (M ~ 1O 13 M and a ~ 100 kpc) 
we find d/d (s 7 ^ m ) ^ 5 • 10~ 2 . Taking an hypotetical very massive and compact object, 
with M ~ 10 13 M Q and a ~ 10 kpc we get d/d (s r ^ m ) ^ 0.2. Thus, a significant neutrino 



conversion effect in the galactic halo is excluded, in contrast with the result in ref. [TJ 



Conversely, significant matter-induced conversion can be realized in halos surrounding 
a cluster of galaxies. Taking the mass of a cluster as M = (10 13 -j- 1O 15 )M and the size 
a ~ 1 Mpc, we obtain from eq. (|92j) 

d < 10 29 + 3 • 10 32 cm" 2 ~ (10- 2 + 20)d {sf n ) . (93) 

For the maximal value , rf/<io(s™ m ) ~ 20, from the condition fl82|) we get the sensitivity 
to the mixing: sin 2 28 ^ 3 • 10 -3 . With this value we find that the adiabaticity condition 
(|1|) is fulfilled for Am 2 ^ 4 • 10" 6 eV 2 . 

The maximal sensitivity is achieved in the energy range of the Z-resonance, where 
also inelastic scattering and absorption are important. Indeed, in section 4.4, taking 
8^26* = 0.3, we have found that d a b s ~ d min at s ~ 6-10 3 GeV 2 , where d min ~ 3-10 32 cm~ 2 . 
This value coincides with the upper edge of the interval ([53|). For smaller sin 29 d m i n is 
larger, and the absorption effect on oscillations becomes even more important. Notice 
that do(sz) takes its minimum value at the Z-resonance: for neutrino energies outside the 
resonance d (sz) is larger, so that d/d (s z ) < 1 and the minimum width condition is not 
satisfied. 

Thus, we have found that in the halos of clusters of galaxies a significant matter- 
induced conversion can be achieved in the narrow interval of energies of the Z-boson 
resonance, where, however, the absorption and the effects of inelastic interactions are 
important. 



5.5 Early Universe 

In this section we consider neutrinos produced by cosmologically distant sources and 
propagating in the universe. The refraction occurs due to the interaction of the neutrinos 
with the particle background of the universe made of neutrinos electrons and nucleons. 
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The number densities of baryons, n&, and electrons, n e , are given by nj = n e = 
TjbTij, where t]b = 10~ 10 4- 10 -9 is the baryon asymmetry of the universe and n 7 is the 
concentration of photons. At present time n 1 = n° ~ 400 cm" 3 . We will describe 
the neutrino background by the total number density n u + n D and the CP-asymmetry 
r] u = {n v — n p )/n T The value of r\ v is unknown. A natural assumption would be r] u ~ f]f,: 
in this —almost CP-symmetric— case the total concentration of neutrinos equals n v + n v = 



4n 7 /ll [39]. However strong asymmetry, r) v ~ 1, is not excluded. The upper bound 



r] u ~ 10 for muon and tau neutrinos follows from the Big Bang Nucleosinthesis |40|,[41 . 

For T] u ~ 1 the contribution of the neutrino background to refraction dominates and 
the interaction of neutrinos with electrons and nucleons can be neglected. It can be 
checked that the contributions of neutrino-electron and neutrino-nucleon scattering to 
the refraction width are smaller than do at any time after the neutrino decoupling epoch, 

tdec — Is. 



In the framework of the standard Big-Bang cosmology, the number density of neutrinos 
in the universe decreases with the increasing time as 



n(t) 



n 



t > t 



eq 



n 



eq 



(94) 



t < t 



eq 



Here to — 10 s is the age of the universe, and t eq ~ 10 s is the time at which the 
energy densities of radiation and of matter in the universe were approximatively equal. 
We denote by no and n eq the neutrino concentrations at t = to and t = t eq respectively. 

The matter width d(t) crossed by the neutrinos from the time t of their production to 
the present one is given by the integration of the concentration (p4|): 



d(t) 



to 



n(r)dr 




t > t, 



eq 



(95) 



t < t 



eq ; 



where dy = t n is the present width of the universe and d eq 

dbti t — tp 



du [to/ t eq — 1] is the width 



"eg- 



In what follows we will focus on the case of matter domination epoch, t > t eq , for 



which the width (|9"5|) can be expressed in terms of the redshift, z = (t /t) 



2/3 



eq- 

as: 



d(z) = d\j (z + 1) 



di 



[z + D 



(96) 



where di = tn = du(z + 1)2 is the width of the universe at the time t of production of 
the neutrino beam; n is the concentration of the neutrino background at t. According 
to eq. (0), for large enough z the width at the production time, d iy gives the dominant 
contribution. 
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Another important feature of the propagation of neutrinos from cosmological sources 
is the redshift of energy. The refraction width do depends on the center of mass energy 
squared s of the incoming and the background neutrinos. As a consequence of the redshift, 
Sz = Sz(z), the width d changes with time (with z) during the neutrino propagation: 
do = d (sz(z)). Thus, the width of matter d should be compared with some effective 
(properly averaged) refraction width do which in fact depends on the channel of transi- 
tion. For non-relativistic background neutrinos of mass m„ we have that s increases with 
z as: s ~ 2m v E oc (1 + z), where E is the energy of the neutrino beam. For relativistic 
background with energy Eb one gets s ~ 2EbE oc (1 + z) 2 . 

Let us consider neutrino propagation in a strongly CP-asymmetric background, 3> 
ni, with r\ v ~ 1. For simplicity we assume also flavour symmetry: n Ve = n v = n„ r . 

For the v a — v s channel the refraction width ( |68| ) increases smoothly from its low energy 
value, d (sz < 1) = 3c/ /4, to the high energy one, d (sz 3> 1) = d . For neutrinos 
produced with energy E ^ 10 21 eV and mass of the background neutrino m Va ~ 2 eV we 
get Sz ^ 0.5, which undergoes redshift during the neutrino propagation. Therefore we 
can use the low energy value of do as the effective refraction width. From fl9"6| ) we get the 
ratio r = d/do(sz) in terms of the redshift z and the asymmetry r\ v : 

r{z) = ~ 2.2- lCrVO + 1) 1 • (97) 

Taking the maximal allowed asymmetry, r] u ~ 10, we find that r = 3 is reached at z ~ 5, 
which corresponds to rather recent epoch. Possible sources of high energy neutrinos, the 
quasars, have been observed at such values of the redshift. With smaller asymmetries, 
r\ v & 1, the minimum width condition requires much earlier epochs of neutrino production, 
z^27. 

Notice that in general the minimum width condition, d > do/ sin 29, is not sufficient to 
ensure a significant transition effect: as discussed in section 2, the width d\i% required to 
have conversion probability P > 1/2 depends on the specific effect involved. For the adi- 
abatic conversion in varying density (section 2.3) we have found the result di/ 2 — 1.5cf min 
(see eq. (^)). Therefore, the condition P > 1/2 requires larger values of r(z), r ^ 4.5. 
From eq. fl97D , with r\ v ~ 10, we get r ~ 4.5 for z ~ 7.5. 

For v a — u s channel the refraction width do, eq. (|69D, has a resonance character with 
absolute minimum do{s 1 ^ m ) ~ do/7 at s^ in ~ 1. Outside the resonance it takes the values 
do{sz < 1) = 3d /4 and d (sz S> 1) = d . For neutrino energy E ^ 10 21 eV at production 
and mass of the background neutrino m Ua ^ 1 eV we get sz < 1, so that we can use the low 
energy value of the refraction width and the result for the ratio r coincides with that in eq. 
(|97|). For neutrinos produced at time t = U with extremely high energies, E ~ 10 21 4- 10 22 
eV and mass of the background neutrino m Va ~ 1 4 3 eV we get sz > 1 at the production 
time, so that, due to redshift, sz(z) will cross the resonance interval, in which d has 
minimum. This, however, does not lead to larger values of the ratio r(z), since the time 
interval At during which sz remains in the resonance range is short: At/tj ~ l-5jz — 0.04. 



29 



The matter width collected in the interval At is d res ~ dAt/ti ~ l.b'jzd, and the ratio 
r(z) for the resonance epoch equals r(z) = d res /do{s^ m ) ~ d/4do- That is even smaller 
than r outside the resonance, eq. (^); therefore the result (|97|) holds also in this case. 

Thus, in the extreme condition of very large v — v asymmetry and production epoch 
at z ^ 5 the matter width crossed by neutrinos satisfies the minimum width condition. 
Let us comment on the character of the matter-induced neutrino conversion in this case. 
After its production, the neutrino beam experiences a monotonically decreasing density. 
Moreover, the energy of neutrinos decreases due to redshift, which also influences the 
mixing. Taking into account the decrease with time of both the neutrino energy and 
concentration, the adiabaticity condition (|18|)-(|19|) can be generalized as: 

lit) =1ot = 70(^ + 1)^ < 1 (98) 
to 

70 ~ 3K,t tan 2 2# ' (99) 

where Vq is the present value of the neutrino-medium potential. For the present epoch we 
get 70 — 10 2 / tan 2 29, so that the adiabaticity is strongly broken. For tan 2 26 ^ 0.1 the 
adiabaticity can be realized at t/t < 10~ 3 , or z > 10 2 . Thus, for z < 10 2 the conversion 
takes a character of oscillations in the production epoch. A detailed study of the dynamics 
of neutrino conversion in the universe will be given elsewhere [JO . 



For the flavour channel v a — vp the refraction effect appears at high energies, sz ~ 1, as 
a consequence of different energies of the background neutrinos v a and vp. The absolute 
minimum of the refraction width (|7T|), d nm {s l z ) = 3do, gives a value of r(z) which is 4 
times smaller than that in eq. fl97|). Correspondingly, even for the extreme conditions 
of r\ v ~ 10 and z = 5 we get r ~ 0.8. That is, significant matter effect is excluded for 
neutrinos from the oldest observed sources. The value r ^ 3 can be reached for z <; 14. 
Notice that, according to eq. (|7T|) , the refraction disappears for low energies, s z 1, 
and the redshift spoils the conditions of absolute minimum of do, even if it is realized in 
certain epoch. Therefore, larger z are required to have significant conversion effect. 

Similar conclusions can be obtained for v a — vp channel, where the refraction width 
([74]) has the local minimum do(s z ) ~ 67^0 due to Z- resonance. This minimum is re- 
alized however during the resonance epoch At. Taking the corresponding matter width 
d res — l-5jzd, we get r(z) = d(z)/4do, which is even smaller than in the v a — vp case. 



Let us consider CP-symmetric neutrino background, rii = rq, with the assumption of 
flavour symmetry: n Ue = n u = n Ur . As discussed in section 4.4, unsuppressed refraction 
effect appears for active-sterile neutrino channels at high energies, s z ~ 1, where the 
propagator corrections become important. 

For v a — v s (and similarly for v a — v s , due to CP-symmetry) the refraction width 
(|76l) has the absolute minimum c? (s™ n ) ~ 2'jzdo, realized at s z im = 1 ± 72, eq. (|77|) . 
Assuming that the neutrinos are produced just before the resonance epoch, we find that 
the width collected during the interval At equals d res ~ dAt/ti — l-S^zd, where d is given 
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by eq. ©, with d v = 2n%/ll 



7 • 10 29 cm 2 . For the ratio r(z) we find: 



r 2 



3d(z) 



2 ■ 10"^ + 1)3 



(100) 



which is significantly smaller than r(z) for u a — u s channel in CP- asymmetric background, 
eq. (]97). With the maximal redshift z ~ 5 eq. (|100|) gives r ~ 0.03, which excludes matter 
effect for neutrinos from the most distant observable sources. The result ( |100| ) holds also 
for the u a — Up (u a — up) channel, since the refraction width (79) has analogous behaviour 
to the one for the u a - 
resonance. 



case, eq. (|76|), with the same local minimum do(s z ) ~ 2^y z do at 



In conclusion, we have found that the matter effect for neutrinos crossing the universe 
is mainly due to the neutrino background. For neutrinos from observable sources (z ~ 5), 
significant conversion effect can be achieved in the u a — v s and v a — v s channels, if the 
background has strong CP-asymmetry, close to the maximum value, r\ v ~ 10. The matter 
effect for the other conversion channels and for the CP-symmetric case is suppressed as a 
consequence the redshift. 



6 Conclusions 

1) . Matter effects can lead to strong flavour transition even for small vacuum mixing 
angle: C 1. This however requires a sufficiently large amount (width) of matter 
crossed by neutrinos: the minimum width condition, d > d min , should be satisfied, where 
dmin = n/(2y2Gp tan 29) = d / tan 29, for low neutrino energies, s <C Mf, and conver- 
sion probability P > 1/2. The absolute minimum d m i n is realized for uniform medium 
with resonance density n r e es . 

2) . We have shown that for all the other realistic situations the required width, dx/ 2 , is 
larger than d m i n . In particular, we have found that ^ 2 /4in = 1 + (1 — 7r/8)7 2 for oscil- 
lations in medium with slowly varying density (7 C 1); d^/im > 1.5 for conversion in 
medium with varying density; d 1 / 2 /d min = tt for castle wall profile. 

3) . We discussed the minimum width condition for high energy neutrinos. For s ^ the 
minimum width d m i n becomes function of s, due to the propagator effect: d m i n = d m i n (s). 
The function d min (s) depends on the channel of interactions: in the case of W (or Z) 
exchange in the s-channel d m i n (s) decreases in the resonance region by a factor ~ 20 with 
respect to the low energy value: d min (0) / d m i n {Ml^) ~ 20. In this region, however, the 
inelastic interactions become important, damping the flavour conversion. 

4) . As a case of special interest we have studied the refraction of high energy neutrinos 
in neutrino background, which can be important for propagation of cosmic neutrinos in 
galaxies and intergalactic space. Again we find that the v a — v a annihilation channel gives 



31 



enhancement of refraction at s ~ M§ , so that d m i n can be ~ 1/27^ ~ 20 times smaller 
than that at low energies. In the case of flavour channels the refraction can appear as the 
result of the difference of masses of the background neutrinos, even if the concentrations 
of the various flavours are equal. 

5) . The minimum width condition allows one to conclude on the relevance of the mat- 
ter effect without knowledge of the density profile, once the width d is known. In some 
astrophysical situations the total width on the way of neutrinos can be estimated rather 
precisely (e.g. by spectroscopical methods) although the density distribution is unknown. 
Significant matter effect in excluded if d < do, or d < d min if the mixing angle is known. 

6) . From practical point of view, a study of the matter effects should start with the check 
of the minimum width condition, d > d m i n . This condition is necessary but not sufficient 
for strong conversion effect. If it is fulfilled the ratio d/d allows one to estimate the mini- 
mal mixing angle for which significant transition is possible: sin 29 > d /d. This condition 
gives an absolute lower bound on 9, which can be achieved for the case of uniform medium 
with resonance density. In other words, given the width d of the medium, the highest 
sensitivity to the mixing angle is achieved if the matter is distributed uniformly and the 
density coincides with the resonance value for a given neutrino energy. For media with 
non-uniform matter distribution the sensitivity to 9 is lower. The stronger the deviation 
from the constant density, the lower the sensitivity. 

7) . We applied the minimum width condition to neutrinos in AGN and GRBs environ- 
ment. For AGN the width d can be estimated by the experimental data on the X-ray 
spectrum, without assumptions on the the density profile. We got d/d Q ^ 1CT 10 for radi- 
ally moving neutrinos, strongly excluding matter effects. In the case of GRBs the width 
d can be evaluated under the assumption that the object is transparent to protons. We 
found d/do ^ 1CT 5 . Therefore, no significant conversion is expected either. 

8) . For neutrinos crossing the halos of galaxies and clusters of galaxies the matter effect 
is given by the interaction of neutrinos with the neutrino component of the halo. We 
have found that for galactic halos the minimum width condition is not satisfied: the re- 
sult d(halo)/do < 0.1 excludes any significant conversion effect. For halos of clusters of 
galaxies we got d(halo)/d Q 10, and the minimum width condition can be satisfied for 
large enough mixing: sin 26* <: 0.1. 

9) . We have considered the refraction of neutrinos from cosmologically distant sources, 
interacting with the neutrino background of the universe. Significant active-sterile conver- 
sion can be expected in case of large v — v asymmetry. We have found that for r\ v = 0(1) 
the condition diuniverse) /d ^ 1 can be achieved for neutrinos from sources, galaxies or 
quasars, with redshift z > 5. The effect on detected neutrino fluxes from these sources 
could be a distortion of the energy spectrum. 
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